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1. Introduction 

In the last years there has been an increasing interest for the con- 
struction of finite-dimensional dynamical systems from soliton equa- 
tions, through the so-called methods of stationary flows and restricted 
flows (see [0], [0 and references therein). The discovery of suitable 
sets of coordinates has allowed one to write the reduced systems as 
physically interesting Hamiltonian systems. In the case of the KdV hi- 
erarchy, the ^-representation for stationary flows has given rise to the 
Henon-Heiles system M, the square eigenf unctions representation 
for restricted flows has furnished the Neumann and the Gamier sys- 
tems |5], . However the relation between dynamical systems which are 
obtained through different reduction techniques from the same soliton 
hierarchy is not clear; moreover a systematic way to find the second 
Hamiltonian formulation for stationary flows of any order, without the 
use of a Miura map, is still lacking. 

The aim of this paper is to give a contribution in these directions. In 
particular: 

i) A bi-Hamiltonian formulation for stationary flows of the KdV 
hierarchy in a suitably extended phase space is derived in a sys- 
tematic way. As an example, the bi-Hamiltonian structure of 
Henon-Heiles-type systems is explicitly shown. 

ii) A map between stationary and restricted flows of the KdV hierar- 
chy is obtained, based on the generating function of the Gelfand- 
Dickey (GD) polynomials. As an application, a map between an 
integrable Henon-Heiles system and the Gamier system with two 
degrees of freedom is constructed. 

iii) An integrability criterion is proposed, which can be applied to 
both stationary and restricted flows. Though weaker than the bi- 
Hamiltonian formulation, it does not require the extension of the 
phase spaces. 

The paper is organized as follows. In Sect. 2 we construct the 
stationary flows associated to the the KdV hierarchy through the kernel 
of the Poisson pencil. Using the generating function of GD polynomials 
as in 0, we give a bi-Lagrangian and a bi-Hamiltonian formulation of 
the Lax-Novikov stationary equations of any order; as an application, 
we exhibit a generalized Henon-Heiles system. 

In Sects. 3-4 we formulate the method of restricted flows in terms 
of the Poisson pencil instead of the spectral problem as in 0, @. This 
formulation allows us to explicitly construct a map between restricted 



ON THE INTEGRABILITY OF STATIONARY AND RESTRICTED FLOWS 3 



and stationary flows, by means of an appropriate extension of the corre- 
sponding phase spaces. The previous map is specialized to the Henon- 
Heiles and the Gamier systems. 

In Sect. 5 we show that the entire bi-Hamiltonian hierarchy of the 
Henon-Heiles and the Gamier systems cannot be reduced from the ex- 
tended to the standard phase space. For this reason, we propose an 
integrability criterion holding for a generic finite-dimensional Hamil- 
tonian system. It generalizes the criterion introduced in M for the 
particular case of the Henon-Heiles system. Though weaker than the 
bi-Hamiltonian scheme, it assures Liouville-integrability of a Hamil- 
tonian system in its standard phase space, i.e. without the intro- 
duction of supplementary coordinates. This criterion is applied to the 
generalized Henon-Heiles system and to the Gamier system with two 
degrees of freedom. 

Now we give some preliminaries, mainly to fix notations and termi- 
nology. Let M be a n-dimensional manifold. At any point u G M, the 
tangent and cotangent spaces are denoted by T U M and T*M, the pair- 
ing between the two spaces by 
<,>: T*M x T U M -> R. For each smooth function / G C°°(M), 
df denotes the differential of /. M is said to be a Poisson manifold if it 
is endowed with a Poisson bracket {, } : C°°(M) x C°°(M) -> C°°(M), 
possibly a degenerate one; the associated Poisson tensor P is defined 
by {/, g}(u) :=< df(u), P u dg{u) >. So, at each point u, P u is a 
linear map P u : T*M — > T U M, skew-symmetric and with vanishing 
Schouten bracket | I0| . A function h G C°°(M) with a non trivial dif- 



ferential df G KerP is called a Casimir of P: P u df(u) = 0. A map 
$ : M — > M is a Poisson morphism if {f,g} o $ = {/ o $,g o $}, 
for each /, g G C°°(M); $ leaves invariant the Poisson tensor P: 
P<s-{u) = &*P U &*, where and $* denote, respectively, the tangent 
and the cotangent maps associated to $. In particular, if the Poisson 
bracket is non degenerate, i.e. if P is invertible, and the Poisson mor- 
phism is a diffeomorphism, $ defines a symplectic (canonical) trans- 
formation. M is said to be a bi-Hamiltonian manifold if it is endowed 
with two Poisson tensors Pq and P\ such that the associated pencil 
P x := P\ — XPq be itself a Poisson tensor for any A G C [ITT], |l2|. 



2. Stationary flows and Henon-Heiles systems 

2.1. KdV hierarchy and Gelfand— Dickey polynomials. Let M 

be a bi-Hamiltonian manifold: if the associated Poisson pencil P x : = 
Pi — XP admits as a Casimir a formal Laurent series h(X) 
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KX) := E hi > (2-1) 

i>o 

then /i is a Casimir of Po an d the coefficients /i, (j > 1) are the Hamil- 
tonian functions of a hierarchy of bi-Hamiltonian vector fields Xf 



Xj = P x dhj = P dh j+1 (j > 0) . (2.2) 

At any point u 6 M, the bi-Hamiltonian flows are given by du/dtj = 
Xj(u), tj being the evolution parameter of the jth flow. The vector 
fields Q2.2p are Hamiltonian also with respect to the Poisson pencil P x . 
In fact the recursion relation ( |2.2| ) can be written as 



Xj = P x dh {j \X) , h u \\) := (X J h(X)) + , (2.3) 

where the index + means the projection of a Laurent series onto the 
purely polynomial part. 

Let M be the algebra of polynomials in u,u x ,u xx , . . . (u = u(x) 
is a C°° function of x and the subscript x means the derivative with 
respect to x), and let Po and Pi be the two Poisson tensors of the KdV 
hierarchy 

„ d „ d 3 d j i 

P :=^, P x := — + Au— + 2u x . 2.4 

dx dx 6 dx 

The gradients of the Casimirs of the associated Poisson pencil P A can 
be obtained searching for the 1-forms v(X) := J2j>o v j which are 
solutions of the following equation 

B\v(X),v(X)) = a(X) , (2.5) 

where a(X) = J2j>-i cijX'i, ctj are constant parameters and B x is the 
bilinear function 

B x (wi, w 2 ) ■= w lxx w 2 + WxWixx - w Vx w 2x + (Au - X)w 1 w 2 . 

(2.6) 

In fact B x is related to the Poisson pencil through the relation 

^-B x (w u w 2 ) = Wi P x w 2 + w 2 P x wi (Vwi, w 2 ) . 

dx ^ 

Eq. ( |2.5| ) can be solved developing the left hand side as a Laurent series 
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B x (v(X),v(X)) = J2 B ^ k (2-8) 

fc>-l 

so that, for each a(A), it furnishes the coefficients of the solution v(X) 
(unique up to a sign). The solution corresponding to a(A) = —A is the 
so-called basis solution u(A); its first coefficients are: 



Vq — 1 , Vi — 2u , v 2 — 2(u xx + 3m 2 ) , f 3 = 2(u ( ^ 4 - ) + 5w 2 + lOu^w + 10m 3 ) 

(2.9) 

and so on, namely the gradients of the first KdV Hamiltonians. In the 
following we shall consider also the 1-form v(X) = c(X)v(X), which is 
solution of (p.5|) for 



a(A) = -Ac 2 (A) , c(A) = 1 + £ c J X ~ j > 

j>i (2.10) 

where the coefficient Cj are free parameters. In this case the first 1- 
forms of the hierarchy are Vq = 1, v i = V\ + ci, v% = v<i + C\V\ + C2 and 
so on. 

The coefficient in ( j2.8p can be expressed through the GD poly- 
nomials. For each Laurent series v(X) let us consider the functions 
£( fc )(A) := B x (v(X),v (k \X)), where «<*>(A) := (A fc w(A)) + ; these func- 
tions have the form 

fl(*)(A) = A fc+ \ 2 + £ A^(p Qi - v v j+1 ) + ]T A~^ (j, k e Nk) • 
3=1 i>o (2.11) 

It can be shown that 



B-i = -Vq , B k = p ok - v v k+1 (k E Nk) ; 

(2.12) 

furthermore, if v (A) is a solution of Eq. ( |2.5|) , the coefficients pjk in 
(|2.11| ) are polynomials in u and its x-derivatives. They will be referred 



to as Gelfand-Dickey (GD) polynomials and the function B x as their 
generating function. 

The fundamental property of the GD polynomials, stemming from 



( |2.11| ), fl2.5|) , (|2.7|) and (|2.3|), is the following relation with the gradients 



Vj = dhj and the bi-Hamiltonian vector fields 
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-T^Pjk = VjX k . (2.13) 

We report some GD polynomials to be used in the following (vq = 1): 
Poo = Au - v i , 

p 01 = Suv 1 -v\-v 2 + 2v lxx , 

P02 = 4wf + 8uv 2 - 2v 1 v 2 -v 3 - v\ x + 2v x v Vxx + 2v 2xx , 

p 12 = 8uv!V 2 -vl + 4uv 3 - ViV 3 -v 4 + 2v lx v 2x + 2v 2 v lxx + 2v 1 v 2xx + \k x 

p kk = 2v kxx v k - v\ x + Auv I . 

The GD polynomials corresponding to the basis solution v(X) are the 
polynomials defined in Prop. 12.1.12]. 

2.2. The method of stationary flows. The method of stationary 



flows [0, [14], [15|] was developed in order to reduce the flows of the 
KdV hierarchy onto the set M n of fixed points of the nth flow X n of 
the hierarchy: 

M n := {u | X n (u, u x , . . . , u {2n+l) ) = 0} . (2.15) 

As M n is odd-dimensional it cannot be a symplectic manifold; nev- 
ertheless we will show that it is a bi-Hamiltonian manifold: it will 
be referred to as extended phase space. Moreover, M n is naturally fo- 
liated, on account of ( |2.2| ) and (|2.4|), by a one-parameter family of 
2n-dimensional submanifolds S n given by 

S n := {u | v n+1 {u,u x ,... ,u {2n) ) = c} (2.16) 

(c being a constant parameter), which are invariant manifolds with 
respect to each vector field of the KdV hierarchy, due to the invariance 
of the 1-forms v k . So M n can be naturally parametrized by v±, . . . , t> n +i 
and by their x-derivatives V\ X) . . . , v nx . We shall use these coordinates 
in the following. 

Here we perform two different stationary reductions of the KdV flows 
by improving the procedure introduced in [0. On one side, we choose 
as a reduction submanifold just the leaf S n of the foliation ( [2.16| ) 



corresponding to c = 0; it is a level set of the GD polynomial po n , 
due to Q, Q and fl2T2|). On account of Eq. (|2TT3| ), also the 



GD polynomials pj n , restricted to M„, are invariant with respect to 
each flow of the hierarchy; thus we can choose as a second reduction 
submanifold a level set of p nn . The one-parameter family of the 
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level sets of p nn forms a foliation of the manifold M n different from the 
previous one. Finally we construct the bi-Hamiltonian structure in the 
ground manifold M n . 

From the computational point of view, one proceeds as follows. 



i) Due to ( p.3| ) and ( |2.5| ), the manifold M n is defined by the solutions 
u of the equation 

B x (v(\),v (n) (\)) = A n a(A) , (2.17) 

where v(X) = J2]=o v j^~^ a W = YjfL-i OjA~ 3 °. In particular if a(A) = 
— Ac 2 (A), as in (|2.10|) , M n is given by 



M n =\u\X n + Y. c jXn-j = \ , (2.18) 

3=1 



i.e. by the solutions of the Lax-Novikov equations |L3|. Taking into 
account Eq. ( j2.11| ) and choosing a_\ = —1, equating in Eq. ( |2.17| ) 
the coefficients of A n+1 we get Vq — 1; from now on we put Vq — 1. 
Moreover equating the coefficients of the other powers of A we get the 
following system: 

Pok - v k+ i = a k (k = 0, . . . , n - 1) , p jn = a n+j (j = 0, . . . , n) . 

(2.19) 

ii) In order to obtain the first Poisson tensor P , we eliminate from 
Eqs. (|2.19| ) u = t>i/2 + a /4 using the first equation (k = 0) and we 
extract the system of n second order ODE's in the Vj (j — 1, . . . , n): 

Pok - v k+ i = a k (k = 1, . . . , n - 1) ; p 0n = a n , 

(2.20) 

which will be referred to as Po~system. The remaining equations ( |2.19| ) 
will furnish a set of n independent integrals of motion. In order to 
obtain a second Poisson structure, we consider the following system 
(Pi-system) 

Pok - Vk+i = a k (k = 1, . . . , n - 1) ; p nn = a 2n , 

(2.21) 

with u as above. 
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iii) The system ( |2.20| ) can be written in Lagrangian form. To this 



purpose, we use the so-called Newton or r-representation introduced 
in []i~6| . Namely, we choose as new coordinates in the first n coef- 



ficients Tj of the formal series r(A) := Jv(X) 



rk = A_, 



(^KA)) (k = l,...,n), (2.22) 

where A*, means the coefficient of X k in a Laurent series. Taking into 
account Eq. (|2.17| ), and observing that 2r n+i = —J21=i r j r n+i-j, Eqs. 
( |2.20| ) are equivalent to 



(A"(r- + (n + ^)r-^)J + = 0. (2.23) 
This system is Lagrangian, with Lagrangian function 

Lf = A_ (n+1) (£(A;V(A)) , (2.24) 
where C(\; ^!(A)) is given, for each Laurent series w(X), by 

£(A; 3(A)) := _( 3| (A)) - ^(3. + — )^(A) - ^ 



The Lagrangian gradients ^ := ^- — -^q^- of are 

^ = ^ I (A-(-r.-(r I + ^)r + ^)) + (* = l,...,^. 

We remark that it is also possible to put also the Pi-system fl2.21|) 
in Lagrangian form. To this purpose, we take as coordinates in 
Qk = r k (k = 1, • • • , n — 1) and g n = \/—v n . By this choice the system 
(|2.21|) is equivalent to 



fez + (?i + -j-)q n - TT = ' ( ^ 2,27 ' ) 
4 4 9„ 



where (A n 1 ?(A)) := (A n 1 y / t>(A)^) + . This is a Lagrangian system 
with Lagrangian 
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L« = A_„( £ (A; U(A)) + f U| - |(U„ + ^ - ^ ^ ^ 
Indeed it can be verified that the Lagrangian gradients of are 

W " A*_! (V^-fa, - (a + ^f*)g + A)l (A; = 2, . . . , n — 1 



;2.29) 



oq k \ ^4 Aq 6 

~5qZ = ~ qnXX ~ + T 5 ** + AqJ 

The two previous Lagrangian systems can be put in canonical Hamil- 
tonian form. For the Po-system the canonical momenta are s n+ i_ re = 
r kx {k = 1, . . . ,n) and the Hamiltonian function 

flf ) = A_ (n+1) (W(A; V(A), /(A))) , (2.30) 
where s(A) = £? =1 SjA -3 ' and 7i(A; 3(A), t(A)) is given by 

W(A; 3(A), ((A)) = ^t 6 (A) + f + ^)3 e (A) + ^ . 

For the Pi-system the canonical momenta are p n = q nx , p n ~k — qkx (k = 
1, . . . , n — 1), and the Hamiltonian function is 



HiP = A_„(H(A; U(A), ^A))) + | ^ + |(n. + -)U^ + ^ 

withp(A) =E"=iPiA-^. 

The two Hamiltonian functions depend, respectively, on the two sets 
of coordinates and momenta (r K , s^), (qk,Pk) and on the two sets of free 
parameters (do, • • • , a n -i) a n) and (ao, • • • , a n -i, «2n)- 

iv) Now let us consider the manifold M n (|2.18| ), which can be parametrized 
either by (r re ,s re ,a n ), or by (qk,Pk,d2n), with a n and a 2 „ as additional 
dynamical variables in M n . On this manifold one can extend trivially 
the canonical Poisson structures, the Hamiltonians and the vector fields 
associated with each one of the two systems as in [IT]. In particular 



the vector fields can be extended in such a way that they are tangent 
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to one of the foliations and Sj$ . Taking into account, on one side, 
the relation between the two sets of coordinates through the original 
variables (v k , Vk x ), on the other side the relation between the two inte- 
grals of motion a n and a2 n through the GD polynomials po n and p nn , 
a map $ : M n -> M n , (r fc , s k , a n ) i-> (q k ,p k , a 2n ) can be systematically 
constructed. It relates the Hamiltonians and the vector fields of one 
system with the corresponding ones of the other system. Since this map 
is not a Poisson morphism, the extended canonical Poisson structures 
associated with one chart is mapped into a Poisson structure different 
from the extended canonical structure associated with the other chart. 
If this second Poisson tensor is compatible with the extended canonical 
one, a bi-Hamiltonian formulation of the two systems is obtained. 
In conclusion we can state the following : 



Proposition 2.1. The Po-system ( P.20T ) and the Pi-system (|2.21|) 



written respectively in the coordinates r k and q k are natural Lagrangian 
systems. The corresponding canonical Hamiltonian systems 



dH^ dHj® 

r kx = ^ ) s kx = o > (2.33) 

uSfe or^ 
opk oq k 



have n integrals of motion given by 



K 3 = ~nVjn, Y = a n+j (j = 1, . . . , n) , Hj = ~-Pjn ]x = a n+j (j = 0, . . . , n - 1) . 
8 8 (2.35) 

Moreover, the map $ : M n — > M n in the extended phase space generates 
a second Poisson structure. 



Remark 2.1. The symbols \Y and \X in (|2.35|) mean that, in the GD 
polynomials Pj k , the coordinates {v k ,v ka ) must be replaced by the 
canonical coordinates (r k ,s k ) and (qk,Pk) respectively and that the 
first order x-derivatives of momenta must be eliminated by means of 
the Hamiltonian dynamical equations (2.33), (2.34). □ 

In the next Subsection we shall give some applications of the results 
stated in this proposition. 

2.3. The bi— Hamiltonian structure of a Henon— Heiles system. 

We consider a generalized Henon-Heiles system with two degrees of 
freedom. Its Hamiltonian is 
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Ha 



UpI+pI 



04 



2 



2 i 1 2 
?1 + |?2 



Oi 

T 91 



(2.36) 



where qi,q2,Pi,P2 are the canonical coordinates and momenta and ao, 
Oi, 04 are free constant parameters. This Hamiltonian encompasses the 
two cases a = a 4 = and a = a\ 
H is related with the Hamiltonian 



introduced in 181. Moreover 



H H=2 (Pi + Pi) + 2 + 5g 2 2 ) + ^ + + g^2 



04 

(2.37) 

through the map 



qi = q [ + 25 , q 2 = q' 2 , a = -2A + 125 , a x = -A 2 + 16AB - 485 2 . 

2 " (2.38) 

The function Hh is the Hamiltonian of a classical integrable Henon- 
Heiles system [|19j with the additional term a\j%q 2 2 . 

The function ( j2.36j ) is the Hamiltonian of the the vector field ob- 
tained reducing X (u) = u x to the stationary manifold M 2 given by 
the fixed points of the flow X 2 + c\Xi + c 2 X n 



M 2 = |u|u (5) + I0u xxx u + 20u xx u x + 30u x u 2 + ci(u xxx + Qu x u) + c 2 u x = o} , 

(2.39) 

where c\ = —a n /2, c 2 = — Oi/2 + Oq/4. 

It can be obtained specializing to the case n = 2 the Hamiltonian 
(|2.32| ) of the Pi-system. In this case H% = Hq and the canonical 



coordinates and momenta are, respectively, q\ = vi/2,q 2 = \/—v 2 , 
Pi = Qix, P2 = <l2x- The integrals of motion obtained by the reduction 
of the GD polynomials are 



u - 1 u - 1 ° 4 

^0 = -gP02, x , ti 2 = -~P22 }X = — g , 

Hi = ~-Pi2 ]x = P 2 qi - PVPi<h - -q x q 2 - -q 2 + ~ -jQiQ 2 + t|#) 

The corresponding Hamiltonian vector fields will be denoted by Xj+x := 
E dHj (j = 0, 1, 2); E being the canonical (4 x 4) Poisson matrix. The 
Henon-Heiles vector field X% is: 
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i 1 2 ^1 ^ OiQ iT 



Xi = [pi,p 2 , -3g x - -g a ~ a o<?i + — > -9i9a + ~ -rQs] ■ 

1<I - 4 (2.41) 



4 



The second Hamiltonian formulation can be obtained specializing to 
the case n = 2 the Hamiltonian ( |2.30| ) of the Po~system: 



rr(O) 5 4 5 2 1 2 a 3 3 2 , a 2 a i 

H y 2 = sis 2 - g^i + 2 r i r 2 - ~~ Y Tl 8 aiTl a ° riT2 ~ ~4 Tl ~ ~ r2 



(4.42) 



where the canonical coordinates (|2.22|) and momenta are, respectively, 
f\ = Vi/2, r 2 = v 2 /2 — vf/4, s\ = r 2x , s 2 = r± x . The integrals of motion 
obtained by the reduction of the GD polynomials are 



K 



Ko = — 



1 

"oP02|y 

■P22, Y 



a 4 , * 



1 



(0) 



-5^2 + sis 2 ri + -s\ - -r\ + 2r x r2 



Oi 3 a 2 ai «o 2 a 2 2 

H ri r,r 2 -\ rir 2 H r 9 — 

4 1 2 1 2 2 2 



a 2 



(2.43) 



and the corresponding Hamiltonian vector fields will be denoted by 



Y; 



EdKj. 



Now we construct the bi-Hamiltonian structure of the Henon-Heiles 
system. Let M 2 be the 5-dimensional extended phase space parametrized 
by (ri, r 2 , s\, s 2 ; a 2 ) or (qi, q 2 ,pi,p 2 ; a^j. It is convenient to make use of 
block notations. So, for example, we denote with (r, s; a) the 5-tuple 
(ri,r 2 , si, s 2 ; a 2 ), with X = [X r , X s ; X a ] T the generic vector field and 
with dK = [dK/dr, dK/ds; dK/da] T the generic gradient of a function 
K (the superscript T means transposition). In this notation a vector 
field X = P dK with Hamiltonian function K with respect to a Poisson 
tensor P will be written 



- X r- 




' prr 


prs 


pra' 


X s 




psr 


pss 


psa 






par 


pas 


paa 





' dK ' 




dr 




dK 




ds_ 




dK 




- da - 



(2.44) 



where P sr = —(P rs ) etc From the definition of r±, r 2 and q±, q 2 in 

terms of V\ and v 2 , and from (|2.40|) and ( |2.43|) one obtains the following 
map $ : M 2 — > M 2 , (r, s; a 2 ) h-> (q,p; a 4 ) 
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Ol 
Pi 



Ol 
P-2 



;-2r 2 -r?) 1/2 
si +ris 2 



32, P2 = -7— ^ ^2TT7^ ' fl 4 = ~ 8K 2 (2.45) 

(— 2r 2 — rf J 1 /^ 

with X 2 given by Eq.( 2.43 ). In these two charts let us consider the 
extended Hamiltonians Hj and Kj, the vector fields Xj (XJ = XJ, 



Xf 



X?, X? = 0) and Yj (YJ = YJ, 



Yf = Yf, Yf = 0), the exten- 
sion of the canonical Poisson structure, E := (-100V The following 

V / 

proposition holds 



Proposition 2.2. The action of the map $ : M 2 —* M2 defined by 
(2.45) on the Hamiltonians Hj, the vector fields Yj and the Poisson 
tensor Pq := E is given by <&*(Hj) = Kj, $*(Yj) = Xj and by 



P :- 





—A T 
8(X!) T 



where A 



J_ ( - 92 \ D _ J_ / -pa \ 

q%\-q2 2<ll J > q%\P2 ) 



A 
B 

8(x?r 
.□ 



-8X|' 
-8Xf 




(2.46) 



Thus we have recovered in the extended phase space M 2 a second 
Poisson tensor Pq. We can check that Pq is compatible with Pi = 
E. Furthermore Pq and Pi give rise to the following bi-Hamiltonian 
hierarchy 



X j+1 := Pi dHj = P dH j+1 {j = 0, 1) , (2.47) 
the Hamiltonians Ho and if 2 being Casimirs of Pq and Pi respectively. 



3. Restricted flows and Garnier systems 



The method of restricted flows was introduced in [20] as a non lin- 
earization of the KdV spectral problem and was generalized in 0. 3 . 
We formulate this method putting the emphasis on the role of the GD 
polynomials and of their generating function; this formulation allows 
us to construct a map between stationary and restricted flows in the 
next section. In view of the applications, we begin by applying the 
method to the KdV hierarchy, recovering the Garnier system. 

Let us consider the following system 
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Poo-^i = «o, Po(^i-E^) =°> P Xk Pk = (fc = l,...,n) 

(3.1) 

where: Ai, . . . , A n are distinct fixed parameters, P Xk := P± — X k Po {Po 
and Pi being the two KdV Poisson tensors (p.4p ). This is a system 
of (n + 2) equations in u, i>i, . . . , (3 n . The second equation will be 
referred to as the P^-restriction of the first KdV flow Xq = PqVi = v± x , 
and the last n equations define the kernel of n Poisson tensors extracted 
from the Poisson pencil. On account of (2.14), ( |2.4| ) and( |2.7| ) this 
system is equivalent to the following one 



«i = Eft+c, B x *({3 k ,[3 k ) = f k , 

i=i (3.2) 



where c and f k are free parameters and B x is just the generating func- 
tion (|2.6| ) of the GD polynomials. 

Using the first two equations to eliminate u and v i from the last n 
equations, one gets a system of n ODE's of second order for (3\, . . . , (3 n : 



2(3 kxx (3 k - (3 2 kx + 2(3 2 k (J2 Pi + d )- = fk (k = 1,. . . ,n) , 



i=i 



(3.3) 



where d := c + ao/2. Introducing the so-called eigenfunction variables 
= J~Pj and the momenta Xj = tyjxi Eqs. (|3.3|) can be written in 
canonical Hamiltonian form 



G 



9Xj 

with Hamiltonian 



ok: 



G 



dipj 



(j = 1,... ,n) 



(3.4) 



1 n 1 

z i=l s 



/ n \ * ra n _p 

va=i / j=i i=i 



The corresponding Hamiltonian vector field 3^g — £ ^ s 



(3.5) 



oo 



oo 



% = [X|, "T 1 ^ + VIM + ^(A| - G + ~7r\~=>] T (I = oo, . . . , \) , 

^ (3.6) 



ON THE INTEGRABILITY OF STATIONARY AND RESTRICTED FLOWS 15 



£ being the (2n x 2n) canonical Poisson matrix. Eqs. ( |3.4| ) are just 
the equations of the Gamier system with n degrees of freedom || . A 
set of integrals of motion is 



X 2 i + 



k=l 



k=i 4 V 



with YTj=\ Ij — 2/Cg- These integrals were obtained in [23] by means 
of a Lax representation; we shall recover them in the next section by 
the use of the generating function of the GD polynomials. 

Let us consider the (2n + 1) extended phase space .M e with coor- 
dinates (ipk,Xk',d) an d the extended Hamiltonian JCq, the vector field 

yn = £ diCn with £ — ( -L o™ o ) . In this space the Gamier system 

V o oo/ 

has a second Hamiltonian structure given by 



(3.7) 



Pi := 



A — ip®ib Ay% 

-(A - ib <g> ip) T x®^ ®X 4^ G 



<y *) T 



<y 6 l 



(3. 



where ® denotes the tensor product, ip — [ipi, . . . , ip n } T , X — [Xi> • • • > Xn] 
A = diag(\i, . . . , A n ). This structure is an extension of the one con- 
structed in for /fc = , (k — 1, . . . , n). In view of the applications we 
specialize the above structure to the case n = 2, in the five-dimensional 
extended phase space M. € with coordinates (/0i, 4>2, Xii Xi'i d)- The fol- 
lowing proposition holds: 



Proposition 3.1. The Gamier vector field 3r 
following bi-Hamiltonian hierarchy 



yc belongs to the 



y j+1 = Vi dQ j = V dg j+1 (j =0,1), (3.9) 
where the Hamiltonians Qj are given by 

Go = ^, & = -(A x + A 2 )^ + i(/i + / 2 ) 

G 2 = AiA 2 ^ - i(A x + A 2 )(/ 1 + h) + i(Ai/! + A 2 I 2 ) , (3.10) 

Q and Q2 being Casimirs of Vo and Vi respectively, and Ii,/ 2 being 
the extensions to A4 e of the integrals of motion ( |3.7| ). 
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As in the case of the Henon-Heiles system, a bi-Hamiltonian struc- 
ture for the Gamier system seems to naturally exist only in its extended 
phase space. Nevertheless in Subsect. |5.3| a realization of the integrabil- 
ity structure introduced in Prop. |5.1| will be constructed in the original 
four-dimensional phase space. 

4. A MAP BETWEEN STATIONARY AND RESTRICTED FLOWS 

Now we shall construct a map between the nth stationary flow and 
the previous restricted flow of the KdV hierarchy. To this end we extend 
the corresponding phase spaces, regarding some free parameters in the 
Hamiltonian functions as additional dynamical variables. 

4.1. The general case. As for the Pi -formulation of the stationary 
flow (2.34) we extend its phase space to a (3n+ l)-dimensional space, 
M n , with coordinates (q k ,Pk] a>o, ■ ■ ■ ,CLn-i,d2n)', analogously we con- 
sider the Po~formulation of the first restricted flow ( |3.4[) in the extended 
space M n with coordinates (if) k , Xk] fi, ■ ■ ■ ,fn,d). 

Let us consider the solutions q k of the dynamical equations (2.34); 
then t>( n )(A) given by 



,(») 



(A) = A g 2 (A) 



(n-1) 



(4.1) 



with q(X) = 1 + Y%=iQj^ K satisfies (|2.17|) and consequently the fol- 
lowing equation 



P A (V n) (A)y n) (A)) = X 2n a(X) , (4.2) 

where, as above, we put u = v\/2 + a /4. So, for each n-tuple of 
distinct complex parameters Xj, any solution t>^ n ^(A) of Eq. ( [4.2|) fulfills 
the system 



B x "(v^(X k ),v^(X k )) =A 2 "a(A fe ) (k 



,n) 



(4.3) 



where v^ n \X k ) 



,(n) 



(A) 



|A=A fc 



In order to have a solution satisfying 



also the second equation (|3.2|) , the Lagrange interpolation formula can 
be used [[Ti], |22fl . It allows us to represent the polynomial v ( n \X) by 



A 



xr ' 



(4.4) 



where p(X) = YYj=i (A — Xj), and 
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/3 fc = ^¥ (k = l,...,n), (4.5) 

( p'(A) means the derivative of p(X) with respect to A). Obviously the 
n functions /3& (|4.5|) are solutions of the following system 



2(3 kx A - PL + 2(3l(± b + %-± A,) - \ k ft = (k = 1, 

j=l 1 j=l WW) ( 4 g) 

furthermore (3 k satisfy the so-called Bargmann constraint 



E (ft - = «i , (4.7) 

3=1 



as one can verify by means of ( |4.4j ). Comparing ( |4.6| ) with ( |3.3j ), we 
can state the following 

Proposition 4.1. Lei ^ : M n .M„, (q,p; a , . . . , a n _i,a 2n ) i-> (ip,x, fi, ■ ■ ■ ,fn,d) 
be the map: 

, _ / £"=<} g=0 gigj-f A fc~ J - gn \V 2 _ Ej=l (Ij^lPn-l -QnPn 



Xk 



(p'(A fe )(E-= 1 ELo^Ar-g 



2n+l 



\P \ A k)j j =0 j=n+l Z j=l 



ao A . (4.< l 



(k = 1,... , n), where Hj are the Hamiltonian functions ( 2.35 ). If 



(qk,Pk) ore solutions of the stationary flows (2.34), then (ipk,Xk) are 
solutions of the Gamier system (|3.4j) for fk and d given by 



Remark 4.1. The function B x is also a generating function of integrals 
of motion for the Gamier system. Indeed evaluating the function B x by 
means of (|4.4j) and eliminating the first ^-derivatives of Xk by means 
of the Hamilton equations ( |3.4j) , one gets 

^iVg(A4r +2d -^S§ (49) 



where Ij are the functions ( |3.7| ). Taking in this equation the residues 
at A = A,- it follows that the functions L are integrals of motion along 



the flow (PIP . □ 
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4.2. The map between the Henon— Heiles and the Gamier sys- 
tem. Now we specialize the map of Prop. 4.1 to the Henon-Heiles 
and the Gamier systems with two degrees of freedom: we obtain the 
surprising result that the Henon-Heiles vector field is mapped into the 
Gamier vector field. Let us consider the seven-dimensional phase space 
of the Henon-Heiles system M 2 with coordinates (q,p; a , a 1; a 4 ). Sim- 
ilarly, for the Gamier systems let us select the parameters fx, f 2 , d and 
enlarge the phase space to a seven-dimensional phase space M. 2 , with 
coordinates (tp, x\ hi hi d). It is easy to prove the following 

Proposition 4.2. Letm : M 2 -»• M 2 , (q,p; Oq, «i, a 4 ) i-> (V>,XI h, h,d) 
be defined by 

^ = + 2A igi - ql) 1 ' 2 , 4> 2 = \ u 1/2 (-\l - 2A 2?1 + qlf' 2 
_ (Aipi - q 2 p 2 ) _ (Aggl ~ ?2P2) 

Xl ~~ , vl/5 I X2 — . >1/9. ) 



(Ai 2 (A? + 2A igi - q 2 2 )) 1 ' 2 ' " (A 12 (-A| - 2A 2gi + ql)) 1 ' 2 

h = Ki (-Xl + a Xi + aiXl - 8H \j - 8H 1 \ 1 + aA , ( 4 - 10 ) 



h = Xn (-X 5 2 + a Xi + ai \l - 8H Xj - 8H 1 \ 2 + a 4 ) , d = ^ - (A x + A 2 ) 
where Ai 2 = Ai — A 2 . The tangent map ^ maps the extended Henon- 



Heiles vector fields Xi, X 2 (|2.47|) into the extended Gamier vector fields 

yi,y 2 



(Xt) =3^, (x 2 ) = y e . (4.ii) 

Moreover the pull-back of the Gamier integrals of motion Goc one? Q e 
are integrals of motion for the Henon-Heiles system 

nG°°) = -l(Xl + Xl) + ^(X 1 + X 2 ) + ^ 
V*(Ge) = Ar 2 2 (2AiA 2 # + (Ai + A 2 )#! + 2H 2 ) + 

+ ((A? + X\) - f (A? + A 2 ) - y(Ai + A 2 )) (4 ' 12) 

The action of the map ^ on the Poisson tensor E of the Henon-Heiles 
system, furnishes a new Poisson tensor for the Gamier system com- 
patible with £. Moreover the action of ^ on the Poisson tensor P is 
given by 
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12 






A 





A r 


B 















where 
A 

B 



^V 2 

v^2 2 



+ + Ax - A 2 



'IplM^l +^2 



-faM^l + ^) V>M + V>2 + A 2 " Ai) 





'(X2^1 - Xl^2 



(X2^1 - Xi^j) 




(4.13) 



(4.14) 



So the map \I/ is not a Poisson morphism. However, according to 
Eqs. ( |4.11| ), the orbits of the Henon-Heiles system are mapped into 
the orbits of the Gamier system. 



5. A NEW INTEGRABILITY STRUCTURE 

5.1. The reduced structures of the Henon-Heiles and the Gar- 
nier systems. In order to have a bi-Hamiltonian hierarchy also in the 
original phase space for the Henon-Heiles and the Gamier systems, one 
can try to apply the reduction techniques known from the literature 



1Q|, |24fl . In particular, two methods can be followed: a restriction to 
the standard phase space or a projection onto it. However, in both 
cases, these attempts fail. 

As for the Henon-Heiles system, if the restriction submanifold is 
chosen to be a leaf S$ of the second natural foliation in M 2 , the 
Hamiltonians Hj, the vector fields Xj and the Poisson structure P\ can 
be trivially restricted respectively to Hj, Xj and E; but it turns out 
that Po cannot be restricted. So two integrable Hamiltonian vector 
fields are obtained in Sj^ but not a bi-Hamiltonian hierarchy. 

If n : M 2 -> S 2 , (qi,q2jPi,P2;a 2n ) >-> (ft, 92, Pi, P2) is the projection 
map, the Hamiltonians Hj and the vector fields Xj cannot be projected 
onto 5*2, because they depend on the fiber coordinate. Instead, the 
Poisson tensors Pq and P\ are projected onto: 



p H := n,p rr 





-A T 



A 
B 



n.pn* = e 



(5.1) 



with A, B as in Prop. 2.2. Because these operators are compatible and 
invertible, one obtains the following Nijenhuis tensor [|25] 
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N H := PjjE = 

and consequently the hierarchy of Poisson tensors Pk '■— N^Ph ,k G Z. 
However these tensors are not invariant along the flow of the Henon- 
Heiles vector field X\ ( |2.41| ). In other words Xi is neither a symmetry 
of P nor of P\ , so that these tensors cannot generate a bi-Hamiltonian 
hierarchy starting from X\. 

As in the case of the Henon-Heiles system, one cannot reduce the 
bi-Hamiltonian structure of the Gamier system with n degrees of free- 
dom onto the restricted phase space. If n : A4\ — > S^°\ x\\\ D l— *■ 

("0H, X\\) 1S the projection map, the Poisson tensor Vo and V\ are pro- 
jected onto two compatible tensors 



A 
B A T 



(5.2) 



A -ip ® ip ' 

-(A - ip ® ip) T x®^ ®X 

(5J 



They give rise to the Nijenhuis tensor Mg := VgS 00 together with the 
hierarchy of Poisson tensor fields V\\ '■— Mg £ , \\ G Z. However these 
tensor fields are not invariant along the flow of the Gamier vector field 
yg ( |3.(j| ), so they do not generate a bi-Hamiltonian hierarchy starting 
from yg. 

5.2. A new integrability criterion. In the previous subsection we 
have put into evidence some problems arising in the geometrical reduc- 
tion of a bi-Hamiltonian structure from an extended phase space onto 
the original one. As an alternative construction, here we introduce 
a new integrability scheme, weaker than the bi-Hamiltonian one, but 
living in the standard phase space. We shall define this new structure 
for a generic Hamiltonian system with n degrees of freedom; for n = 2 
it coincides with the one introduced in |8j for the Henon-Heiles system 
with the Hamiltonian ( 2.37]) and = 0. As new examples of this in- 



3) 



tegrability structure, the case of the Gamier system with two degrees 
of freedom will be discussed here whereas multidimensional extensions 
of the Henon-Heiles system will be presented elsewhere. 

Proposition 5.1. Let M be a In- dimensional Poisson manifold equipped 
with a Poisson tensor Qq, and Zq a Hamiltonian vector field with 
Hamiltonian ho: Zq = Qo dh$. Let there exist a tensor M : TM — > TM 
and a skew-symmetric tensor Q\ : T*M — > TM such that 

Qi = MQo . (5.4) 
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Denote by Zi and cti the vector fields and the l-forms obtained, re- 
spectively, by the iterated action of the tensor M on Zq and its adjoint 
TV* : T*M -> T*M on a := dh 

Z i: =Af%, ai :=M* l a Q (i = l,... ,n-l) . 

(5.5) 

Let there exist n — 1 independent functions hi (i — 1, . . . , n — 1) and 
(n 2 + n — 2)/2 functions (i = 1, . . . , n — 1; < j < i) with fi 0Q = 1, 
Ha 7^ (z = 1, . . . ,n — 1), such that the l-forms oij can 6e written as 

i 

a i = ^2fi ij dh j (i = 1, . . . ,n - 1) . (5.6) 

j=o 

Under the previous assumptions the following results hold: 

i) t/ie vector fields Z{ satisfy the recursion relations 

Z l+ i = Q a i+1 = Qxcti (i = 0,...,n-2). (5.7) 

ii) the functions hi are in involution with respect to the Poisson bracket 
defined by Qq and they are constants of motion for the fields Z^ 

{/ii,MQo = 0, £*„«>)=/, (5.8) 

where Cz t denotes the Lie derivative with respect to the vector field Z/, . 

iii) the Hamiltonian system corresponding to the vector field Zq is 
Liouville-integrable. In addition if Q\ is a Poisson tensor field, then 
also Z\ is an integrable Hamiltonian vector field and the functions hi 
are in involution also with respect to the Poisson bracket defined by Q x . 

Proof. 

i) From Eq. (|5.4j ) and the skew-symmetry of Qo and Qi it follows that 
QoN* = MQ and QiW* = MQi. Then 

Z x - g «i = Z X - Q Af*a = Z X - TVQoao = (5.9) 
and the first relation (|5.7|) is proved by induction since it is 



Zi+i ~ Qo"i+i = - QqM*^ = Af(Zi - Q ai) 
The second relation ( |5.7| ) follows from 



(5.10) 
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Z i+X - Qm = MZi - Q x on = M{Z t - Q oai ) . (5.11) 

ii) By (|5.6|) , the gradients dhk can be expressed for any k in terms of 
dho 

k 

dh k = CEvkiN^dho , (5.12) 

where v\.% are the elements of the matrix a -1 , a being the lower tri- 
angular matrix defined by = (i > j), a^- = (i < j), = 
0, . . . , n — 1). Thus 

{hi,hj} Qo : = (dhi,Q dhj) 

1 j b 
= ViaVjb(N* a dh , Q M* dh ) 



a=0b=0 (5.13) 

i 3 



= T,T,^ jb (dh 0l Af a+b Q dh ) 

a=0 b=0 

and the first relation ( |5.8D follows from the skew-symmetry of the ten- 
sor Af m Qo for any m. Furthermore 

^|(()) = (dhi,Q Q a k ^i) 

k 

= (dh^Qo^fMkjdhj) 

> ** < 514 > 

= zZ^kj{hh hj} Qo 
= 

iii) Since Zq is a Hamiltonian vector field, it is Liouville-integrable on 
account of the previous result. Moreover, since it is 

{hi,hj} Ql : = (dh^Qxdhj) 

i i b 
= Y,Y,^ j b^* a dh ,Q 1 Af* dh ) 

a=0b=0 (515) 

i 3 



J2T,^ jb (dh ,Af a+b Q l dh ) 
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it follows that if Qi is also a Poisson tensor, {, }g 1 is a Poisson bracket, 
Z\ is a Hamiltonian vector field and then it is Liouville-integrable. □ 

Remark 5.1. The recursion scheme and the integrability of the vector 
field Zq do not require that the skew-symmetric tensor Qi be a Poisson 
tensor; so M is a Poisson manifold, not a bi-Hamiltonian one. □ 

In view of the applications of the next subsection, it may be worth- 
while to remark that the results of Prop. |5.1| hold true if the role of Qq 
and Qi are interchanged; to be more precise, one can prove (just as for 



Prop. [53 



Proposition 5.2. The integrability scheme of Prop. is still valid 
if Qo is skew-symmetric, Qi is a Poisson tensor and the role of Z Q 
is now played by Z\ = Qidh . The involution relations ( p.8| ) become 
{h i ,h j } Ql = 0. 

5.3. The integrability structure of the Henon— Heiles and the 
Garnier systems. In Subsect. |5J] we have recovered by projection 
onto the quotient manifold S 2 the Nijenhuis tensor ( |5.2|) and a hierarchy 
of compatible Poisson tensors; however, it is not possible to associate 
to these tensors and to the Henon-Heiles vector field X\ ( |2.41| ) a bi- 
Hamiltonian hierarchy of vector fields. Nevertheless it is possible to use 
these elements to construct an example of the integrability structure 
introduced in Prop. 5.2. To this purpose, let us make the following 
choices: 

i) Qi = E, the vector field Z\ := X\ Q2.41D with Hamiltonian ho '■= 
H ( gg§ ); 

ii) the tensor field Af : 
P H as in Q; 

iii) the function h\ := H% ( 2.4U|) and the functions /i^- as /iio 



N, 



H 



and Qc 



P-5 



Nx P H , with 



0,/iu = \jq%\ 

then it is immediate to check that the conditions of Prop. 5.2 are 
satisfied. Moreover the vector field Z := Q dh = P_ 2 dif is a new 
integrable vector field: 



-2pxqx - p 2 q 2 
-PiQ2 



-p 2 2 + 6gr + 2gig 2 
P1P2 



CL4 

4q* 



Zaoqf + f q\ 



Y + ?>q\q2 - — q2 + aoqii2 



(5.16) 



This integrability structure is related, through the map ( |2.38| ), to the 
one introduced in M for the Hamiltonian ( |2 . 3 7| ) with = 0. 
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For the Gamier system with two degrees of freedom one can con- 
struct an example of the integrability structures of Prop. |5.1| . Indeed 
if one uses the elements of Subsect. 5.1 and makes the following choices: 

i) Q := £, h := Q x (^10|) , Z := y g (pj) ; 

ii) N := Mq 1 = 8Vg°° } with Vg as in (5.3) , := V-oo = M g °° 8 



A 1 A 2 



iii) the functions h x := Q 2 ( |3.iq) , /i 10 = 0, /j u = -a^+Ax^-AiAs' 
then the conditions of Prop. |5.1| are satisfied. Moreover the vector field 



:= /-tn^e is a new integrable vector field ( 3^e is the restriction to 
the submanifold of -M g , d = cost, of the vector field 3^2 ( |3.9| )). 

At last, we compute how the map between the standard phase spaces 
of the Henon-Heiles and of the Gamier systems, induced by the map 
( |4.10| ), acts on the recursion operators of the previous integrability 



structures. 

Proposition 5.3. Let us consider the map ^ : (qi,q2,Pi,P2) l— * ("01 > "023 Xi-> X2) 

= Ar 2 1/2 (A? + 2A igi - ql) 1 ' 2 , iJ 2 = (X 12 r 1/2 (-Xl - 2A 2?1 + ql) 1 ' 2 

(A1P1 - q 2 p 2 ) (A 2 pi - q2V2) 

Xi - ~ 7TF> ■> X2 



A 12 (A? + 2A igi - q 2 2 )f 12 ' " (A 12 (-A| - 2A 2gi + g 2 2 ^ V2 '(5-17) 

The map ^ relates the recursion operators of the Henon-Heiles and of 
the Gamier systems: ^*Nh = A/^" 00 ^/*. 



6. Concluding remarks 

In this paper we have derived a bi-Hamiltonian formulation for sta- 
tionary flows, and for the first restricted flows of the KdV hierarchy. 
Our approach amounts to searching the kernel of, respectively, the Pois- 
son pencil and n-Poisson structures extracted from the Poisson pencil 
of the KdV hierarchy. In this approach the generating function of the 
GD polynomials plays a relevant role. Moreover it allows us to con- 
struct a map between stationary flows and restricted flows; in the case 
of the fifth-order stationary KdV equation, this map relates solutions 
of the Henon-Heiles system with solutions of the Gamier system. How- 
ever, to obtain these results one must extend the phase space of the 
reduced flows by means of some free parameters naturally contained in 
the corresponding Hamiltonian functions. This difficulty can be over- 
come, at least if one analyzes complete integrability of a Hamiltonian 
system without requiring an explicit knowledge of a bi-Hamiltonian 
structure. To this purpose, we have introduced a new integrability 



ON THE INTEGRABILITY OF STATIONARY AND RESTRICTED FLOWS 25 

scheme in the standard phase space, which implies Liouville integrabil- 
ity of the reduced Hamiltonian systems. For brevity we have applied 
this scheme only to the Henon-Heiles and the Gamier systems with 
two degrees of freedom. Other examples such as Henon-Heiles type 
systems with three and four degrees of freedom, constructed by means 
of the reduction method of Sect. 2, will be discussed elsewhere. 
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suggestions. 
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